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Herbertsmithite (ZnCu3(OH)6Cl2) is often discussed as the best realization of the highly frus- 
trated antiferromagnetic kagome lattice known so far. We employ density functional theory calcu- 
lations to determine eight exchange coupling constants of the underlying Heisenberg Hamiltonian. 
We find the nearest neighbour coupling Ji to exceed all other couplings by far. However, next- 
nearest neighbour kagome layer couplings of 0.019Ji and interlayer couplings of up to — 0.035Ji 
slightly modify the perfect antiferromagnetic kagome Hamiltonian. Interestingly, the largest inter- 
layer coupling is ferromagnetic even without Cu impurities in the Zn layer. In addition, we validate 
our DFT approach by applying it to kapellasite, a polymorph of herbertsmithite which is known 
experimentally to exhibit competing exchange interactions. 

PACS numbers: 71.20.-b,75.10. Jm,75.10.Kt,75.30.Et 



Quantum spin liquids have fascinated physicists for 
decades as this exotic ground state constitutes a novel 
state of matter^. The magnetic moments in a spin liq- 
uid do not order even at extremely low temperature due 
to a high degree of frustration in the magnetic system. 
Typical examples for lattices that lead to frustration of 
antiferromagnetic interactions are triangular, pyrochlore 
and kagome lattices. While experimental realizations of 
quantum spin liquids have long been scarce, in particular 
the discovery^ of the perfect kagome lattice realization 
in herbertsmithite has led to considerable excitement!^. 
In the eight years since the discovery of the S = 1 /i 
kagome antiferromagnet nature of ZnCu3(OH)6Cl2, nu- 
merous experiments have been performed to ascertain 
the spin liquid groun d sta te of herbertsmithite, its prop- 
erties and excitations^^. In particular, measurements 
of the magnetic susceptibility 6 show antiferromagnetic 
couplings of the order J « 17 meV (~ 190 K) and no 
magnetic ordering down to 50 mK. Muon spin rotation 
measurements^ confirm the absence of magneti c or der- 
ing and inelastic neutron scattering experiments^! find 
that fractionalized quantum excitations are present in 
ZnCu3(OH) 6 Cl2. More recently, non-ideality of the re- 
alization of the kagome Heisenberg Hamiltonian in her- 
bertsmithite due to additional interactions and in the 
form of site disorder has been the focus of many stud- 
ies^. While defects within the kagome layer are detected 
in nuclear magnetic resonance^ but not in recent x-ray 
scattering measurements^, Cu impurities on interlayer 
Zn sites seem to play a role^. Low temperature devi- 
ations between theory for the kagome antiferromagnet 
and experimental susceptibilities as well as anisotropics 
in thermodynamic quantities^ p oint t o a small nonzero 
Dzyaloshinskii-Moriya interactio n 1 13 * 14 !. Evidence of this 
interaction h as be en found in electron spin resonance 
measurements^^. However, the experimental and the- 
oretical discussion about the Hamiltonian correctly de- 
scribing herbertsmithite is far from settled. 

Therefore, we undertake an effort to determine the pa- 



rameters of the underlying Heisenberg Hamiltonian using 
all-electron density functional theory methods. We will 
show in this Letter that the exchange coupling constants 
from first principles corroborate that ZnCu3(OH)gCl2 
is a near perfect realization of a kagome antiferromag- 
net with a dominant coupling of J\ — 182 K. However, 
there are small corrections to this picture: A next-nearest 
neighbour coupling in the kagome layer of 0.019 Ji and in 
particular some interplanar couplings between — 0.035Ji 
and 0.029 Ji could actually be relevant for the nature and 
excitations of the spin liquid ground state in herbert- 
smithite. 

We perform density functional theory calculations with 
the full potential local orbital (FPLO) basis setP^ us- 
ing generalized gradient approximation (GGA|^ and 
GGA+U functionals. The exchange couplings, J i; are 
obtained from total energy calculations for different Cu 
spin configurations in supercells of various sizes^. Be- 
fore proceeding to herbertsmithite, we test our methods 
on kapellasite, a polymorph of herbertsmithite, which has 
been investigated before, both theoretically^ and exper- 
imentally^. We use the structure of kapellasite as given 
in Ref. 1221 and determine the hydrogen position by relax- 
ation^. The structure is shown in Figure [l] We create 
two different supercells: A 2 x 1 x 2 supercell with P 1 
symmetry and 10 inequivalent Cu positions with the pur- 
pose of resolving four interlayer couplings, and a 2 x 2 x 1 
supercell with P 1 symmetry that provides symmetry in- 
equivalent 3rd nearest neighbours in the kagome plane. 
Using GGA+U with U = 6 eV and J = 1 eV, we find 
two significant couplings, J\ — — 14.2 K and Jd = 24.0 K 
(see Table [jj. Other couplings like J2 are significantly 
smaller (around 1 K), and we find interlayer couplings 
to be negligible 24 . This is in very good agreement with 
the observation of Ref. 21 that experimental data are 
compatible with a J\-Jd model with Jj = —15.0(4) K 
and Jd = 12.7(3) K. Note that the ratio between J\ 
and Jd depends on the choice of U as in particular the 
antiferromagnetic Jd is inversely proportional to D^l 
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Even more recently, the high temperature series expan- 
sion method was refined to fit both magnetic susceptibil- 
ity and specific heat data, yielding the set of parameters 
Ji = -12 K, J 2 = -4 K, J d = 15.6 KP. Thus, we can 
proceed with some confidence to analyze the Heisenberg 
Hamiltonian parameters of herbertsmithite. 




FIG. 1: (Color online) (a) Crystal structure of kapellasite, 
viewed along the c direction, (b) Kagome lattice formed by 
the Cu sites in (a). Note that Jd and J4 exchange paths both 
correspond to a distance of 6.3 A but while J4 points precisely 
along a nearest neighbour bond, Jd cuts diagonally across a 
Cu hexagon with nonmagnetic Zn in the center. 



TABLE I: Exchange coupling constants for ZnCu3(OH)eCl2 
(kapellasite) determined from total energies of five different 
spin configurations in a 2 x 2 x 1 supercell. 



name 


dcu-Cu 


type 


Ji (K) 








U = 6eV 


Ji 


3.15 


kagome nn 


-14.2 


Ji 


5.45596 


kagome 2nd nn 


-0.7 


Ji 


6.3 


kagome 3rd nn 


-0.3 


Jd 


6.3 


kagome 3rd nn 


24.0 



We use the structure of ZnCu3(OH)6Cl2 (herbert- 
smithite) with R 3m space group determined by Shores 
et aiP which is shown in Figure [2] (a). A big difference 
with respect to the polymorph kapellasite is that Zn is 
now between kagome layers rather than in the centers of 
its hexagons. In Figure [3j we present the bandstructure 
and density of states. At the Fermi level, we find Cu 
3d states which hybridize with O 2p and CI 3p states. 



As expected, Zn plays no role at Ep. Note that the 
Dirac-point-like feature at K for an energy of 0.2 eV 
above the Fermi level becomes an avoided crossing with 
a tiny gap in a fully relativistic calculation. Based on 
our experience with azurite, another complex quantum 
spin system containing Cu 2+ ions^J and the fact that 
kapellasite, the polymorph of herbertsmithite briefly an- 
alyzed above, was shown to have longer ranged compet- 
ing interactions^, we determine all exchange constants 
up to Cu-Cu distances of 8.6 A. In order to allow for 
determination of the diagonal coupling in the kagome 
lattice, J5 , we double the unit cell along a and prepare a 
structure with Pm space group and 12 inequivalent Cu 
sites. As appropriate for Cu, we employ a GGA+U ex- 
change correlation functional with U = 6 eV, J = 1 eV 
and atomic-limit double-counting correctiorP^l. Total en- 
ergies for nine different spin configurations allow us to 
calculate the eight exchange coupling constants listed in 
Table [TTJ. The three couplings within the kagome layer 
are shown in Figure [2] (b), and the geometry of the three 
most important interlayer couplings is presented in Fig- 
ure [2] (c)-(e). To avoid confusion, we number the cou- 
pling constants Ji strictly according to ascending Cu-Cu 
distances. Within the kagome layer, the most important 
correction to the presently discussed Hamiltonian for her- 
bertsmithite is the next-nearest-neighbour coupling J 3 . 
Theoretical investigations of the J\- J2-kagome lattice in- 
dicate that the nature of the spin liquid g round state 
could depend on such a next nearest coupling 26 2 ' . Mcs- 
sio et aW^ have even extended the range of the couplings 
to the 3rd nearest neighbour across a hexagon; our set of 
parameters would put herbertsmithite in the q = spin 
liquid phase, in agreement with Ref. |9] Very recently, the 
J\- J2-Heisenberg model on the kagome lattice has been 
studied with a pseudofermion functional renormalization 
group method^, very good agreement with the inelas- 
tic neutron scattering experiment of Ref. [9] is reached for 
Ji = 0.017Ji which is very close to our value. 

Now we come to the interlayer couplings. First of all, 
it is important to note that while each Cu site has four 
interactions via J\, four via J3 and six via J5, the in- 
terlayer bonds are numerous: There are four Ji bonds, 
six J4 bonds, eight Jq bonds, six J7 bonds and six Jg 
bonds. Interestingly, we find Ji to be an antiferromag- 
netic interlayer coupling of size 0.029Ji, and J§ a fer- 
romagnetic interlayer coupling of size — 0.035 Ji. Previ- 
ous studies based on the spin-rotation-invariant Green's 
function method showed that a stacked kagome system 
remains short-range ordered independent of the sign and 
strength of the interlayer coupling 2 ^!. We have performed 
a first test of the relevance of the interlayer couplings for 
susceptibility and specific heat using high temperature 
series expansion-^. This method has been very useful to 
discuss the kagome lattice Heisenberg modeP^ as well as 
various additional terms^l. We find that at least in the 
region of applicability of this method the effect of inter- 
layer couplings is noticeable. We hope that our results 
inspire more precise manybody calculations that could 
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establish the consequences of interlayer couplings for the 
low temperature properties of herbertsmithite. 

TABLE II: Exchange coupling constants for ZnCu3(OH)6Cl2 
(herbertsmithite) determined from total energies of nine dif- 
ferent spin configurations. Energies were calculated with 
GGA+U functional at U = 6 eV, J = 1 eV and with atomic 
limit double counting correction. 



name 


dcu~Cu 


type 


J^ (K) 








U = 6eV 


kagome layer couplings 


Ji 


3.4171 


kagome nn 


182.4 


Jz 


5.91859 


kagome 2nd nn 


3.4 


■h 


6.8342 


kagome 3rd nn 


-0.4 


interlayer couplings 


J 2 


5.07638 


interlayer 1st nn 


5.3 


Ji 


6.11933 


interlayer 2nd nn 


-1.5 


Je 


7.00876 


interlayer 3rd nn 


-6.4 


Jr 


8.51328 


interlayer 4th nn 


3.0 


J 9 


9.17347 


interlayer 6th nn 


2.5 


In summary, our ah 


initio-b&sed analysis of the Cu- 


Cu exchanj 


$e coupling 


constants in kapellasite and her- 



bertsmithite provides a detailed description of these ma- 
terials. Our results for the dominant interactions are 
in excellent agreement with experiments. Moreover, we 
are able to resolve the strength and sign of weaker, 
but not negligible, exchange interactions that were not 
known up to now and are important for understand- 
ing the behavior of these materials at low temperatures. 
Both polymorphs, even though they are realizations of 
a perfect kagome lattice, show a few remarkable differ- 
ences. The nearest neighbor Cu-Cu exchange interaction 
is strongly antiferromagnetic in herbertsmithite (~ 190 



K) and weakly ferromagnetic in kapellasite (~ —13 K) 
due to the fact that the Cu-O-Cu angle in herbertsmithite 
is 119° compared to 106° in kapellasite. Kapellasite 
shows a significant antiferromagnetic 4th nearest neigh- 
bor coupling along the diagonal of the Cu hexagon ( Jd) 
which is negligible in herbertsmithite since the exchange 
path in kapellasite is through the in-plane Zn situated 
in the center of the hexagons. Also, the stacking of the 
kagome layers in both polymorphs is crucial for under- 
standing the interlayer exchange couplings. In kapella- 
site, the kagome layers are stacked in a similar fashion as 
in the layered TiOCP^or Cs2CuClJ^ where interactions 
are mostly of van der Waals nature. In this situation, the 
interlayer couplings are comparatively smalP^. In con- 
trast, in herbertsmithite the interlayer Cu-Cu couplings 
are partly through Zn orbitals. This leads to relatively 
significant antiferromagnetic (J2) and ferromagnetic (Jq) 
interlayer couplings. Nevertheless, the ratio between the 
dominant intralayer coupling J\ and the dominant inter- 
layer coupling remains large enough for this system to be 
considered a very good realization of a two-dimensional 
kagome lattice and only at low temperatures should the 
smaller Ji become important. This and the importance 
of couplings other than the dominant ones for the spin- 
liquid behavior in these materials should be investigated 
in the future. 
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FIG. 2: (Color online) (a) Crystal structure of herbertsmithite ZnCu3(OH)eCl2, viewed along the c direction, (b) Kagome 
lattice formed by the Cu sites in (a). Exchange paths between nearest, next nearest and third nearest neighbours within the 
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FIG. 3: (Color online) Band structure and density of states 
of ZnCu3(OH)6Cl2 calculated with GGA exchange correlation 
functional. High symmetry points of the P3m space group 
are M = (1/2, 0, 0), K = (1/3, 1/3, 0) and A = (0, 0, 1/2) in units 
of the reciprocal lattice vectors. DOS is given in states per 
eV per unit cell (containing 3 formula units). 
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SUPPLEMENTARY INFORMATION 



Details for exchange constants of kapellasite 



In Tables SI and SII we provide the results of total en- 
ergy calculations with GGA+U functional using different 
values of U. The 2x1x2 supercell used in the calcu- 



lation for Tabic SII allows resolution of four interlayer 
couplings of kapellasite. They are all very small which 
is not surprising considering the van der Waals gap be- 
tween the layers of kapellasite (see Figure SI). This is 



a significant difference to the polymorph herbertsmithite 
that has kagome layers coupled in the third dimension 
via O-Zn-0 bonds (see Figure S2). 



TABLE SI: Exchange coupling constants for kapellasite de- 
termined from total energies of five different spin configura- 
tions in a 2 x 2 x 1 supercell. 



name 


dcu-Cu 


type 


Ji (K) 
jj — Q e y 


Ji (K) 
jj — j e y 


Ji (K) 
[/ = 8eV 


Ji 


3.15 


kagome nn 


-14.2 


-13.4 


-12.6 


Jz 


5.45596 


kagome 2nd nn 


-0.7 


-0.7 


-0.6 


Ja 


6.3 


kagome 3rd nn 


-0.3 


-0.3 


-0.3 


Jd 


6.3 


kagome 3rd nn 


24.0 


19.8 


16.3 



TABLE SII: Exchange coupling constants for kapellasite de- 
termined from total energies of eight different spin config- 
urations in a 2 x 1 x 2 supercell. Note that here, kagome 
3rd nearest neigbours are symmetry equivalent so that J4 = 
2/3J4 + yaJd. 



name 


dcu-Cu 


type Ji (K) Ji (K) 


J, (K) 






U = 6eV U=7eV £/ = 8eV 


kagome layer couplings 


Ji 


3.15 


kagome nn —11.8 —11.3 


-10.9 


Jz 


5.45596 


kagome 2nd nn —1.2 —1.0 


-0.9 


Jk 


6.3 


kagome 3rd nn 10.8 8.8 


7.3 


interlayer couplings 


J 3 


5.733 


interlayer 1st nn —0.4 —0.3 


-0.2 


Js 


6.54139 


interlayer 2nd nn 0.04 0.04 


0.04 


J 6 


7.91421 


interlayer 3rd nn 0.04 0.04 


0.03 


Js 


8.51806 


interlayer 4th nn 0.01 < 0.01 


< 0.01 



FIG. SI: (Color online) (a) Side view of the crystal structure 
of kapellasite (approximately along b direction), (b) Inter- 
layer exchange paths for the Cu sites in (a). 




FIG. S2: (Color online) Side view of the crystal structure of 
herbertsmithite (approximately along b direction). 



